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INTRODUCTION 


Lunar atmosphere research has tended to center on gases with 
predictably large sources and cm those which have been identified by 

40 

Apollo experiments. An early candidate atmospheric constituent was Ar, 
which was noted by Heyman and Yanlv (1970) to have a surface correlated 
component in returned soil samples, and an abundance in excess of what 
can be explained by potassium decay. The source of the excess argon 
was attributed to atmospheric argon ions which have been accelerated by 
solar wind fields and implanted in soil grains (Manka and Michel, 1970, 

1971) . 

40 

Prior to Apollo 11, Ar was not expected to be an important gas 
on the moon. Most thoughtful predictions included volcanic gases such as 
CO 2 , ^0, and H 2 S, as well as a solar wind supplied component consisting 
mainly of neon. The post Apollo view is that lunar volcanism is essentially 
nonexistent, that the dominant gases of lunar origin are radiogenic helium, 
argon and radon (Hodges, 1977a), and that solar wind helium is the major 
source of lunar atmosphere, albeit most of this helium is trapped in satel- 
lite orbits (Hodges, 1978). Other solar wind elements, notably hydrogen 
and carbon, must exist in the lunar atmosphere, while nitrogen, neon and 
argon from the sun are less certain constituents. 

Hydrogen is an enigma because the large solar wind influx of 
protons is certainly escaping from the moon at essentially the inflow rate, 
despite the fact that no hydrogen has been identified by experimental methods. 
The Apollo 17 orbital UV spectrometer provided an upper bound of 10 hydrogen 
atoms/cc (Fastie et al., 1973) which is grossly inadequate to sustain the 
hydrogen escape rate. Oxides and hydrides of solar wind carbon are found 
in lunar soil samples, suggesting molecular formation within soil grains. 
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and supporting the hypothesis that may be the form of escaping hydrogen. 

Owing to the large influx of solar wind carbon, and the lack of 
an adequate mechanism to continually mix the soil, transporting assimilated 
carbon downward, it has been proposed (Hodges, 1976) that carbon inflow 
and escape are in balance, similar to hydrogen and helium. A major problem 
witu this hypothesis is the lack of experimental evidence for carbon compounds 
in the atmosphere. However, all pertinent data from the Apollo 17 mass spec- 
trometer were obtained during lunar night, when these gases are presumably 
adsorbed on surface materials. There is tentative evidence that CH^, CO, and 
C0 2 exhibit argon-like presunrise increases due to lateral flow of these gases 
being desorbed near the terminator (Hoffman and Hodges, 1975). 

While the distributions and abundances of various carbon gases 
in the lunar atmosphere have not been established by measurement, there is 
good reason to expect that they are present, and that carbon gases dominate 
the bound part of daytime lunar atmosphere. All other candidate gases for 
which model atmospheres have been computed indicate daytime concentrations 

4 

less than 10 /cc (Hodges et al., 1974). 

Early predictions of a large neon atmospheric component biased 
both experiments and data analyses. Data from the Apollo 16 orbital 
mass spectrometer (Hodges et al., 1972) and the Apollo 17 lunar surface 
mass spectrometer (Hodges et al. , 1973) provide upper bounds on nighttime 
neon concentrations which, perhaps fortuitiously, agree with a model atmos- 
phere (Hodges et al. , 1974). Since no good lower bounds exist, the total 
absence of neon at night can noc be ruled out. Its virtual absence from 
the atmosphere would imply only that soil grains are not saturated with 
solar wind neon; this possibility is difficult to refute based on data from 
returned soil sample analyses. In lunar daytime the suprathermal ion 
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detectors (SIDE) at several Apollo sites have recorded a number of events in 
which the solar wind v x B field was briefly aligned to accelerate atmos- 
pheric photoions into one of the SIDE instruments. Generally these events 
are compatible with a neon atmosphere significantly more dense than the 
aforementioned models (cf. Benson and Freeman, 1976). It remains to be 
determined whether a carbon compound, such as CO, or perhaps light metallic 
ions created as the result of solar wind sputtering of the lunar surface 
could also have produced these events. 

The only gases on the moon with realistic, positive identifi- 

40 

cations and long term data bases are helium and Ar. It is important that 
maximum use be made of the available data, both to understand the exospheric 
behavior of argon and helium, and to understand the conspicuous absence of 
several other candidate lunar exosphere constituents from the data. An 
interesting and useful aspect of the argon and helium data it the fact that 
argon is adsorbed at night while helium is not. In addition helium escapes 
kineticly while argon does not. These extreme characteristics assure that 
there is an experimental basis against which to test new theories regarding 
lunar atmosphere dynamics and the fates of the absent species. 

Recent research efforts have been directed mainly toward the 
synthesis of the argon adsorption and diffusion properties of lunar soil, 
and the use of argon-40 as an archetype in attempting to understand the 
migratory behavior of 1^0, mercury halide compounds, and other condensible 
volatiles. Subsequent discussion outlines the most recent progress, 
and includes some preliminary material that is now being prepared for 
publication. 


3 


EXO SPHERE- REGOLITH ENCOUNTER SIMULATION 


A significant result of this program of research was the 

determination by Monte Carlo simulation that following impact on the 

regolith surface of the moon (or Mercury), free exospheric atoms migrate 

vertically in the sol^. in a random manner that is statistically similar 

to a one dimensional random walk. The similarity is that the probability 

of reemergence of an atom from the regolith after a collisions with soil 

-3/2 

grains is asymptotically proportional to n as n becomes large (Hodges, 
1980b). Long sequences of collisions are to be expected in the regolith, 
as are long sequences of "bad luck" in coin flipping games. This is 
necessary to account for the continual downward diffusion of some atoms, 
but it also makes the mean value of n divergent. Hence, the concept of 
a mean interaction time for encounters of exospheric atoms with the 
regolith is meaningless. 

In Monte Carlo simulation of volatile migration on the moon 
an atom is traced through a succession of random thermal-speed ballistic 
trajectories. When a ballistic atom impacts the regolith surface it adsorbs 
on a soil grain, where it resides for a random, temperature dependent time 
interval. Following desorption the atom either emerges from the regolith or 
it collides with another grain where it repeats the adsorption-desorption- 
fLight cycle. Collisions within the soil allow vertical wandering of the 
atom and, as in a one dimensional random walk, the deeper into the soil an 
atom wanders, the more collisions it must make before emergence. The 
probability that an atom emerges from the soil after n collisions is roughly 
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(from an emperical fit of the data in Figure 1 of Hodges, 1980b). In the 
Monte Carlo lunar exosphere simulation, a random number n at each atom- 


(1 v '-^v N 


regolith encounter Is obtained from the inverse of expression (1) 


n - int{1.28/u } (2) 

where the operator ’ lnt* Indicates truncation of the argument to its integer 
value, and u is a random number chosen uniformly from the range 0 to 1. 

Assuming that the atom in question makes its initial contact with 
a soil grain at time t Q , its probability of desorbing at time t^ Is 




where t is the temperature dependent mean desorption time. A random deviate 
of t^ is found by equating expression 3 and a random number u^ chosen uniformly 
between 0 and 1. Neglecting atom flight times between collisions with soil 
grains, the time interval for n adsorption-desorption events is determined 


by the relation 


\ ln(u.) * -ln( n u.) 

■'t T i-1 i-1 


where t is the time of emergence of the atom from the regolith. (Xd-ng to 
n 

the fact that the mean value of the right hand term of (4) is n, with standard 
1/2 

deviation of n , it is convenient to approximate that term by n when n is 


large. 


Integration of the leftmost term of equation (4) requires specifi- 


cation of the temperature dependence of t, and of the time dependence of lunar 
surface temperature. The former is ammenable to synthesis, as wiJ! be dis- 
cussed later, but the latter is somewhat subjective. 
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In previous models of the exospheres of the moon and Mercury it 
has been assumed that radiative equilibrium on a smooth, spherical surface 
is adequate to specify synodic temperature variations (cf. Hodges, 1973, 

1974; Curtis and Hartle, 1978). However, for condensible volatiles the 
strong dependence of desorption time on temperature makes the distribution 
of smaJ 1 scale thermal irregularities too important to ignore. Some high 
lunar surface features are illuminated, and hence heated, up to 100 km 
on the night side of the mean terminator. At the equator this corresponds 
to as much as a six hour advance of local sunrise or delay of sunset. 

Owing to orographic surface slopes that range up to 45°, shadows are 
present at local solar zenith angles greater than 45°. These shadows 
can delay sunrise and adwice cvnset up to 3.7 days at the equator, and 
produce continuous shade at latitudes above 45°. An equatorward tilt of 
a surface increases the solar heat influx. For example, the north wall of 
a typical fresh bowl shaped crater with walls that slope 45° (Wood and 
Anderson, 1978), located at a latitude of 45°N, has the same insolation as 
a level surface a( the equator. The south wall of the same crater receives 
no primary radiation, but it is heated by a substantial infrared flux from 
surrounding sunlit surfaces (Hodges, 1980b). In addition, inhomogenieties 
in subsurface soil conductivity cause differences in nighttime cooling rates, 
producing small scale temperature differences (Schultz and Mendell, 1978) 
that are important in view of the exponential nature of the dependence of 
desorption time on temperature. 

The ideal computer simulator of exospheric volatile migration 
would be based on a thermal model that gives the actual synodic oscillation 
of temperature at any point on the lunar surface. This is impractical, and 


fortunately, unnecessary. What is absolutely needed is a procedure 
to generate an appropriate random set of synodic temperature oscillations 
at any latitude. The scheme devised to do this includes generators to 
represent random distributions of surface slopes, of shadows and sunlight 
in daytime, of the Illumination of tall features on the near-night side 
of the terminators, and of soil conductivity. These parameters are combined 
in emperical formulae to produce synodic temperature oscillations. 

The data on which the surface temperature synodic oscillation 
generator is based were obtained from a large number of calculations of 
the temperature oscillations at various latitudes and for a variety of 
snadow configurations and surface parameters. The method of soil temper- 
ature calculation outlined by Keihm and Langseth (1973) was used, but 
with the modifications that solar illumination was truncated by shadows, 
and a daytime Infrared reradiation flux of 12% of the local Insolation was 
adopted from the crater temperature calculations of Hodges (1980b) to 
account for infrared heating of shaded areas in daytime. This resulted 
in a set of emperical formulae for temperature based on the variable 
parameters of the calculations. In the exosphere simulator, latitude is 
fixed by impact point, while the impact velocity vector direction sets 
bounds on the local slope of the soil. Longitudinal extents of insolation 
are determined randomly from shadow distributions based on the work of 
Watson et al. (1961) and Arnold (1979), and a surface slope distribution 
with a 6° mean slope and 45° maximum. The nighttime cooling rate is 
chcsen randomly from a distribution that causes the model to fit the 
Apollo 17 orbital measurements of small scale nighttime surface temperature 
irregularities reported by Schultz and Mendell (1978). The complete dts- 
crlption of the computation scheme and the emperical formulae used in the 
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exosphere simulator are found in the report entitled 'A Lunar Surface 
Temperature Model for Monte Carlo Exosphere Simulation. 

A key problem with the new lunar surface temperature model is 
that it fails to represent seasonal effects in shaded areas near the 
poles. In the model, latitude is measured with respect to an apparent 
axis that is actually the intersection of the plane of the geometric 
terminator and the plane of the subsolar selenographic meridian. The angle 
between the apparent and true polar axes varies sinusoidally with annual 
period and amplitude equal to the tilt of the lunar axis, or 1.53°. At 
low latitudes this difference is unimportant, but at 85° it causes a + 20% 
variation in the lengths of day and night, and above 87° apparent latitude 
the model is only reasonable at the equinoxes. Owing to the small fraction 
of lunar surface area involved it is convenient to give the polar regions 
ad hoc treatment. What has proved adequate is to ignore differences 
between the apparent and true lunar axes at latitudes below 87®, and to 
assign a trapping time distribution to the polar cap regions^. 

Cold traps present another enigma. In Hodges (1980b) it was 
shown that the^e are permanently shaded areas that are cold enough to 
retain adsorbed argon for long periods, provided the soil grain surfaces 
in the traps are pristine. The work of Arnold (1979) and of Hodges (1981c) 
suggest that water contamination of the traps is unavoidable, making long 
term retention of argon unlikely. An important question that remains 
unanswered is whether the summer polar cap area is complete?'* cleansed of 
adsorbed wat^r and other volatiles, so that in winter it ca*. provide 
temporary traps for argon. 
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INTERACTION OF ARGON ATOMS WITH LUNAR REGOLITH 


One of the most impressive Indicators of the pristine nature 

of lunar soil is the adsorption of exospheric argon-40 at night. Figure 1 
40 

shows the Ar concentration as a function of time through two lunations 

at the Apollo 17 site. The data are roughly restricted to sunset to 

sunrise intervals, while the daytime extrapolations are res nits of an 

exospheric simulation calculation. It can be noted that just past sunset 

the gas concentrations begin to decrease, reaching a minimi m shortly 

before sunrise. This decrease is caused by an excess of adsorption over 

desorption as the lunar surface cools. In the absence of adsorption 

the concentration would have increased at night, varying roughly as the 

-5/2 power of surface temperature (Hodges and Johnson, 1968). Sudden 

desorption of argon at sunrise causes a flux of gas back into night, 

producing the rapid increase in concentration just before sunrise. 

Another important feature of Figure 1 is the distinct decrease 
40 

in the amount of Ar on the moon between March and July, 1973. This 

decrease requires a time varying rate of supply of argon, a phenomenon 

that seems out of character for the otherwise quiescent moon. Speculation 

on the mechanism for sporadic argon releases has ranged from a deep source 

in a semi -molten, core-free asthenosphere of primative lunar material 

(Hodges and Hoffman, 1974 and 1975; Hodges, 1977a) to near surface 

mylonization of crustal rock by thrust faulting (Binder, 1980), and thence 

to sudden desorption of argon from polar cap cold traps due to seasonal 

changes in shadowing (Hodges, 1980b) . 

Recent efforts have been directed toward an improved understanding 

40 

of the interaction of exospheric Ar with lunar regolith. Initial analyses 
of the Apollo 17 mass spectrometer measurements of argon resulted in 
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synthesized adsorption and % minim im sticking tine, which are now 
recognized as mathematical artifices that approximately describe the 
multiple adsorption random walk of exospheric atoms Into and back out 
of the soil (Hodges, 1980b). 

As the argon-regollth process has become better understood, 
the estimate of the activation energy for adsorbed argon has increased 
from 3000 cal/mole (Hodges, 1977a) to 6000 cal/mole (Hodges, 1980b), and 
subsequent discussion will shew that the true value may be as large as 
9000 cal/mole. The activation energy for argon on glass is about 3800 
cal/mole (de Boer, 1968) . The high activation energy on the moon seems 
to be a verification of an intuitively attractive hypothesis, that the 
degree of cleanliness of soil grains on the lunar surface is unattainable 
in the laboratory, and that laboratory activation energy data necessarily 
depict properties of surfaces contaminated with water vapor and possibly 
other volatiles. This Idea seems to be supported by the extremely low 
argon activation energy (500-700 cal/mole) deduced by Frisillo et al. 

(1974) from measurements of BET "C" constants for returned lunar soils 
reported by Holmes et al. (1973). However, the BET "C" constants usually 
give low heats of adsorption because they are determined by the last part 
of monolayer formation, which occurs at sites where activation energies 
are low (Holmes et al. , 1973, and Frisillo et al., 1974). Thus the low 
heats of adsorption of argon deduced from "C" constants may n«., be analogous 
to the adsorption of argon on clean rocks on the lunar surface. The 
difference between laboratory results for argon adsorption on glasr* at _ 
on returned lunar samples Is more puzzling. 

The numerical value of the activation energy is greatly lnfluneced 
by the underlying assumption of functional dependence of desorption time 
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on temperature. Frlsillo et al, [1974] suggested that descrplton tins 
should be represented by 

t - r E/RT (5) 

o e 

where E is the activation energy and t q is a constant having an expected 
-13 

value of 1.6 x 10 s. This expression was adopted by Hodges (1980b), 

resulting in the estimate that E ■ 6000 cal/mole for argon on lunar regolith. 

With the implementation of the more n llstlc surface temperature 

model discussed above, it became apparent that for a constant value of t q , 

equation (5) could not represent adsorption phenomena throughout the lunar 

night. The quantum mechanical theory cf vapor pressure (cf. Kennard, 1938) 

-2 

suggests that at low temperatures t q must vary as T , aid hence that 


_C E/RT 
i 2 


( 6 ) 


is a candidate representation of desorption time. The important difference 
between equations (3) and (6) in the argon exosphere simulator is that (6) 
decreases more rapidly than (5) as temperature increases. 

Figure 2 Illustrates the cype of data now emerging from the 
exosphere simulator calculations. The abscissa i& longitude from the sub- 
solar meridian, but t can also be Interpreted as sun hour angle at a fixed 
point on the rotating moon. In the lower frame the histogram shows the 
computed average flux of argon-40 to the regolith surface at latitudes less 
than 30 s . A line graph superimposed on the nighttime p rt of the plot 
represents smoothed flux data from the Apollo 17 mass spectrometei . The 
upper frame shows the average surface concentration of adsorbed argon, 
obtained by integration of the difference between downcoming and upgoing 
fluxes at the regolith surface. Following sunset, the delayed rise in 
adsorbed argon reflects the slow cooling of the surface, but eventually 
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Figure 2. The downward flux of exospheric argon and the surface concen- 
tration of adsorbed argon a3 functions of longitude measured from the 
subsolar meridian. 


a do min ance of advection over ballistic transport causes the adsorbed 
argon concentration to become nearly constant. The sharp peaks at the 
sunrise terminator, in both the flux and the adsorbed concentration data, 
result from the very high probability that atoms released from the soil 
at the terminator will, in the course of subsequent ballistic flights, 
recycle to the night side of the terminator and be readsorbed. 

In the simulation that produced Figure 2 the cold traps were 
assumed to be contaminated by water vapor and/or other frozen volatiles, 
making argon desorption very fast, and the traps ineffective. However, 
other modeJs, with varying trap conditions, have resulted in similar 
reproductions of the mass spectrometer flux measurements, and in each case 
the surface concentration of adsorbed argon approximately duplicated the 
data shown in the upper part of Figure 2. 

40 

For the rapid presunrise increase in the Ar flux to be reproduced 
by the simulator it is necessary that at the average presunrise surface 
temperature of 92K the desorption time be approximately 36 seconds. This 

requirement fixes the relationship of the parameters C and E in equation (6). 

-17 2 

For the exosphere simulation of Figure 2, the values C * 2x10 sK and 
E * 9306 cal/mole were used. Models with larger values of C tend to give 
lower fluxes following sunset, but a change of C by several orders of mag- 
nitude is necessary to change the post-sunset flux by a factor of two. 

An important descrepancy between model and experimental results 
occurs at the post-midnight minimum of the flux. Here the model gives 
essentially zero because, in the course of a succession of nighttime impact 
events, - it is highly probable that a regolith encounter will occur wherein 
equation (2) generates a very large deviate for the number of adsorptions. 
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When this happens the atom Is retained on the surface until sunrise. 

The near void in the model flux just after midnight indicates that 
almost all atoms suffer this fate. 

The reason for the post-midnight difference between model 
and experimental flux data is probably due to oversis lification in the 
model of the interaction process for atoms that make many adsorption 
collisions. To avoid rapid return to the exosphere these atoms must 
penetrate to depths equivalent to a large number of mean grain diameters. 
At night the very poor conductivity of the upper layer of soil produces 
a large temperature gradient (Keihm and Langseth, 1973), substantially 
decreasing the desorption time for exospheric atoms that penetrate to a 
depth of the order of one millimeter, corresponding to several hundred 
free path lengths. To account for this effect in an approximate way, 
equation (2) has been replaced by 



-4 2/3 

where K is a small number of the order of 10 . For values of u>>K 

the statistical behavior of n given by equations (2) and (7) are essentially 
identical, but with decreasing u the rate of growth of equation (7) becomes 
much slower than for equation (2). Because of the asymptotic treatment of 
large values of n in evaluating the right hand side of equation (4), the 
effect of using equation (7) to generate random values of n is equivalent 
to decreasing the mean desorption time for atoms that penetrate deeply into 
the soil. 

When equation (7) is used to generate random deviates of the 
effective number of collisions of atoms with soil grains the expected result 
occurs: the argon flux throughout the nighttime decay is generally increased 
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due to desorption of atoms ti>at previously would have been retained 
until sunrise. To fit the model to the mass spectrometer measurements it 
is obviously necessary to decrease E and increase C from their previously 
cited values that were appropriate when equation (2) was used to determine 
n. This is gratifying because a much larger value of C is also determined 
from the quantum mechanical theory of vapor pressures, i.e. 


h 3 - 5/2 -12 

c “ w? = 2 ' 3 x 10 


( 8 ) 


where h is Planck's constant, k is Boltzmann's constant, m is atomic mass, 
and o is the effective area of an adsorbed atom. Figure 3 shows model flux 

— 4 

and surface concentration as functions of solar longitude for K m 10 , 

-12 2 

E = 7208 cal/mole, and C * 2 x 10 s K . It must be noted that this 

result is based on an artificial treatment of the temperature gradient in 

the soil, and that there is no present physical justification of the form 

of equation (7). This is an important subject for further research because 

it is pertinent to all phases of volatile migration on the moon. Questions 

40 

raised by the temporal changes in exospheric Ar abundance on the moon, as 
illustrated by the factor of two difference in the two lunations of data 
shown in Figure 1, include whether the source of argon is episodic, and 
how a rapid decay could be compatible with any recycling of implanted argon 
ions (Hodges, 1977b). Pristine polar cap regions that serve as seasonal 

cold traps were proposed by Hodges (1980b) as possible temporary sinks for 
argon, allowing rapid decay of the exospheric abundance following a release, 
and providing sudden releases of adsorbed gas from trap areas that reemerge 
into sunlight with the approach of sumner infrared reradiation levels. 

Ihe data of Figure 4 represent the time histories of two 
40 

impulsive releases of Ar under differing conditions. In part A the winter 


polar caps were presumed to be water saturated and hence ineffective 



Figure 3. Argon flux and surface concentration with soil temperature 
gradient compensation. 
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adsorbers of argon. Permanently shaded cold traps were assumed to 
be ice covered and Ineffective In both cases. The sources of releases 
A and B were at the equatorial subsolar point. In each section of 
Figure 4 the diamond symbols represent the fraction of the released gas 
remaining on the moon as a function of time, the arithmetic difference 
of unity and each diamond being the fraction that has ionized. The X's 
denote the fraction of the release that is in the exosphere and the difference 
bet-reen diamonds and X's is the fraction that is adsorbed on soil grain 
surfaces . 

Primary interest in Figure 4 centers on the exospheric fraction 
data (X's) because these relate to the Apollo 17 mass spectrometer measure- 
ments of argon. In case A the exospheric decay has three distinct phases, 
the first being a rapid decrease from unity at the time of the release to 
about 0.2 in less than 10 days due to adsorption on the night side. The 
second phase appears to last about 100 days and has a time constant of 33 
days. This decay is produced by poleward migration of argon and initial 
capture by the winter polar cap. From 100 days on the exospheric fraction 
and the toted argon extinction have roughly the same time constant of 
2500 days. In case B the second decay phase of the exospheric fraction 
is missing and the third phase has a shorter time constant of 212 days. 

The mass spectrometer measurements of argon-40 on the moon in- 
dicate that the main decay phase of the exospheric fraction should last, 
about 100 days and have a time constant that is less than 100 days. The 
active polar cap model, case A in Figure 4, satisfies this criterion 
absolutely* but the inactive polar cap model does not. In other simulations 
(not shown) where gas was released at high latitudes the second decay phase 
was also missing whether or not cold traps were active. It is not totally 
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Figure 4. Time history of argon-40 released impulsively from the 
equatorial subsolar point. Diamonds give the fraction of the release 
remaining on the moon, while the X's show the fraction in the exo- 
sphere. 
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clear how this relates to the mass spectrometer data, but these preliminary 
results seem to Indicate that some forms of argon cold traps exist on the 
moon as suggested by Hodges (1980b), and that some large argon releases 
occurred at low latitudes during the period of operation of the Apollo 17 
mass spectrometer in 1973. 

What is missing in the exosphere decay data of Figure 4 is 
spatial resolution. Figure 5 shows preliminary data from an exosphere 
model that simulates the time histories of the argon fluxes into 18 
separate mass spectrometers equally spaced on the equator of the moon. 
Inadequate computation time is apparent in the noise that pervades this data 
set. However, the results are adequate to show the initial flow of argon 
into nighttime and the contrasting character of initial sunrises at 
locations that were in day and night hemispheres at the time of the gas 
release. Full implementation of this program is awaiting resolution of 
the outstanding questions discussed above regarding the soil migration and 
desorption processes. 
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Figure 5. Argon flux at 18 equispaeed points on the lunar equator following an 
impulsive gas release at the subsolar point at t=0. 
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